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ABSTRACT 



CO 

o 

tyQ ' We present the analysis of a total of 177 h of high-quality optical time-series photomc- 

Q ! try of the helium atmosphere pulsating white dwarf (DBV) EC 20058-5234. The bulk 

of the observations (135 h) were obtained during a WET campaign (xcovl5) in July 
1997 that featured coordinated observing from 4 southern observatory sites over an 
^> , 8-day period. The remaining data (42 h) were obtained in June 2004 at Mt John Ob- 

■ servatory in NZ over a one-week observing period. This work significantly extends the 

discovery observations of this low-amplitude (few percent) pulsator by increasing the 
number of detected frequencies from 8 to 18, and employs a simulation procedure to 
confirm the reality of these frequencies to a high level of significance (1 in 1000). The 
nature of the observed pulsation spectrum precludes identification of unique pulsation 
mode properties using any clearly discernable trends. However, we have used a global 
modelling procedure employing genetic algorithm techniques to identify the n, I values 
of 8 pulsation modes, and thereby obtain asteroseismic measurements of several model 
parameters, including the stellar mass (0.55 Mq) and T e ff (~ 28 200 K). These val- 
ues are consistent with those derived from published spectral fitting: T c g ~ 28 400 K 
and logg ~ 7.86. We also present persuasive evidence from apparent rotational mode 
splitting for two of the modes that indicates this compact object is a relatively rapid 
rotator with a period of 2 h. In direct analogy with the corresponding properties of the 
hydrogen (DAV) atmosphere pulsators, the stable low-amplitude pulsation behaviour 
of EC 20058 is entirely consistent with its inferred effective temperature, which in- 
dicates it is close to the blue edge of the DBV instability strip. Arguably, our most 
significant result from this work is the clear demonstration that EC 20058 is a very 
stable pulsator with several dominant pulsation modes that can be monitored for their 
long term stability. 

Key words: stars: individual: EC 20058— 5238 - white dwarfs - stars: oscillations - 
~ „ . „ stars: interiors - techniques: photometric 
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1 INTRODUCTION 

Nearly 99% of all stars are predicted to end their lives as 
slowly cooling white dwarfs; their study, among other things, 
provides an important perspective on the active lives of all 
the progenitor objects. In particular, if we can determine 
their internal chemical compositions, we have access to key 
data concerning the products of the nuclear reactions that 
power stars in the previous evolutionary stages. In common 
with other astronomical objects, we are largely limited to 
studying the atmospheres of white dwarfs as this is where 
the photons we detect are created. Thus, detailed spectro- 
scopic observations allow us to measure effective tempera- 
tures, surface gravity values (and therefore stellar masses), 
and atmospheric chemical compositions. From these studies 
over many years we know that the great majority of whit e 
dwarfs divide into two classes (e.g. iMcCook fc SiorJ ll999T ) . 
The DA class have hydrogen atmospheres, and the DB class 
have pure helium atmospheres (no detectable sign of hydro- 
gen and other elements). The DAs account for about 86% 
of all white dwarfs, while the DBs dominate the rest. QQ 

In both of these classes there are objects that provide 
further information that is subtly encoded in stellar flux 
variations. These are the pulsators. The cooler hydrogen at- 
mosphere pulsators (DAV class) were di scovered serendip- 
itously in the late 1960s dLandoltJ 1 19681 ). while their hot- 
ter helium atmosphere cousins were detected following the- 
oretical prediction s and a targeted search (|Wingetl Il98ll: 
IWinget et al]|l982h . 

The use of pulsating stars to infer some of their intrinsic 
properties is called stellar seismology or simply asteroseis- 
mology. Asteroseismology uses the detected pulsation modes 
of a star to constrain computer models and thereby mea- 
sure stellar properties, including otherwise hidden interior 
physical quantities . This endeavo ur has been a very produc- 
tive exercise (e.g. IWingetJ Il998l ) for white dwarfs, in part 
due to the relative simplicity of the white dwarf structure, 
but also due to the potentially rich pulsation mode spec- 
trum that results from the nonradial spheroidal (g-mode) 
pulsation mechanism. All observed white dwarf pulsations 
are attributed to buoyancy-driven g-modes, as the observed 
periods (70 s - 1500 s) are inconsistent with the predicted 
periods of pressure (acoustic) waves in the extremely dense 
white dwarf material: such pulsations would have periods 
several orders of magnitude smaller, and have never been 
observed inspite of a number of search es (|Robinsonlll984 
iKawaler et~al1ll994l ; ISilvotti et ai1l2007l ). 

The more identified pulsation modes, the more success- 
ful is the white dwarf asteroseismology: in effect, each de- 
tected pulsation mode adds another constraint to the mod- 
elling process. 

Given that the DA objects form a substantial major- 
ity of all white dwarfs, it is natural that the DAVs domi- 
nate the known objects for these compact pulsators. In fact, 
following a recent la rge increase in the d etected number 
of faint white dwarfs (|Kleinman et alj [2004l . and follow-up 
programmes to detect more pulsators in these new objects 
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mosphere pulsators has mushroomed to more than 140. 

The helium atmosphere DBV pulsators are far fewer: 
in 2007 there are 17 known obje c ts in this class 
dBeauchamp et"aT. 1 1 19991 : lHandlerl l200ll ; iNitta et al. 1 120051 , 
l2007t) . The southern object EC 20058-5234 (QU Tel, hence- 
forth simply EC 20058) was the 8th DBV to be discovered. 
It has a B magnitude of approximatley 15 and was discov- 
ered by the Edinb u rgh-Cape (EC) fa int blue object survey 
jStobie et al.lll997l ). lKoen et al.1 l|l995l ) first reported a study 
of its properties: they spectroscopically established its DB 
classification and also demonstrated its variability using the 
techniques of time-series photometry. Their 20 hours of pho- 
tometry over a four month period in 1994 revealed a total 
of 8 pulsation frequencies and showed that the object was 
a low amplitude DBV variable that appeared to be quite 
stable; this is in complete contrast (on both counts) to the 
class prototype, CD 358. In fact CD 358 is a large amplitude 
variable which also exhibits considerable changes in its ob- 
served pulsation sp ectrum over variou s timescales ranging 
from days to years (IKepler et al.1120031 ). 

The Fourier analysis of the discovery data set suggested 
the presence of low level frequencies below the chosen sig- 
nificance level, so one of us (P OD) proposed a W hole Earth 
Telescope (WET) campaign (jNather et al.lll990l ) to probe 
for coher ent periodicit ies at lower amplitudes than those de- 
tected bv lKoen et al.l Consequently, a WET run was sched- 
uled for dark time in July 1997, corresponding to the op- 
timal observing season for the target. Following the WET 
run, regular single-site monitoring has been carried out at 
primarly Mt John Observatory (NZ). 



2 OBSERVATIONS 

2.1 WET 1997 photometry 

Observing time covering a period of 9 days on 4 southern 
telescopes was obtained in July 1997 in order to monitor 
EC 20058. Time-series aperture photometry was carried out 
using photometers equipped with blue sensitive photomulti- 
plier tubes and no filters in the light beam at all 4 observing 
sites. Consequently, the resulting 'white light' passbands em- 
ployed in the observations had effective wavelengths similar 
to that of Johnson B, but with a significantly wider pass- 
band. 

All photometers used 10 s integration times and op- 
erated with essentially 100% duty cycles, given the na- 
ture of the photon counting systems employed. The 1.0- 
m telescope at the South African Astronomical Observa- 
tory (SAAO) used a single channel photometer in combi- 
nation with autoguiding; both the 1.0-m telescope at Mt 
John University Observatory (MJUO) in New Zealand and 
the 1.6-m Itajuba telescope in Brazil were equipped with 
two-channel photometers; and a th ree-channel photometer 
l|Kleinman. Nather fc Phillips! fl996l) was used in combina- 
tion with the CTIO 1.5-m telescope in Chile. A journal of 
all the observations is provided in Table Q] 

Good weather at both SAAO and Mt John, combined 
with poor conditions at CTIO, meant that the bulk of the 
observations were obtained at just two observatories. How- 
ever, inspite of the participation of only 4 telescope sites 
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Table 1. Journal of observations of time-series photometry of EC 20058— 5234 obtained during the Whole Earth Telescope (WET) 
extended coverage campaign (XCOV15) in July 1997. Column 2 provides the observatory and telescope for each run, which are: CTIO 
= Cerro Tololo Interamerican Observatory, La Serena, Chile (observer: A. Nitta); MJUO = Mount John University Observatory, Lake 
Tekapo, New Zealand (observer: D. J. Sullivan); OPD = Observatorio Pico dos Dias, Itajuba, Brazil (observer A. Kanaan); and SAAO = 
South African Astronomical Observatory, Sutherland, South Africa (observers: D. Kilkenny, F. van Wyk). Columns 3 and 4 provide the 
UT starting date and time for each run, column 5 gives the total number of 10 s integrations in each run, column 6 gives the run length 
in hours, while the last column is the barycentric Julian start date (BJD - = BJD — 2450000, see text) for each run corresponding to the 
UT start times. Note that the two OPD runs marked with an asterisk in column 1 (ra411 &; ra413) were not used in the final analysis 
due to an unexplained timing error (see text) 
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mation of a pulsation mode that is close to marginal in the 
WET data set, so it is germane to include some analysis of 
these data here. 

The Mt John 1.0-m tel escope in combination with 
a thr ee-channel photometer (jKleinman. Nather fc Phillips! 
1996) using unfiltered 'white light' (and hence a passband 
similar to that mentioned previously) was used to acquire 
time-series photometry of EC 20058 in June 2004. This pho- 
tometer includes the functionality of the one used at CTIO 
in the 1997 WET run, which enables the observer to mon- 
itor 3 regions of the sky using nominally identical minia- 
ture blue sensitive (Hamamatsu R647-04) photomultiplier 
tube s. But, it also incorporates a very useful improvement 
(e.g. ISullivanl 120001 1 . In addition to continously monitoring 
the target star, a comparison star and sky background, the 
comparison star can also function as a continuous guide star. 



overall, the WET campaign achieved an observing duty cy- 
cle of 64% for the whole 8.3 day observing period, and a 
duty cycle of 89% for the central 4.4 days when all tele- 
scopes contributed. The WET campaign complete coverage 
can be gleaned from Fig. [T] where the reduced photometry 
has been graphed in successive 24 hour segments. 



2.2 Mt John 2004 photometry 

Since the 1997 WET campaign, EC 20058 has been regu- 
larly monitored from Mt John during its extended south- 
ern observing season by one of us (DJS), and some 
preliminary reports on th i s work have a lready appeared 
llSuIlivan fc Sulliva 3 l2000l : ISullivanl 120031 ). Quality time- 
series photometry obtained in June 2004 provided confir- 



© 0000 RAS, MNRAS 000, 000-000 



D. J. Sullivan et al. 



12 



15 



18 



UT [hours] 
21 



6 



9 



12 



[2 July 97] 




[3 July 97] 



men 



mm i mm*m$ 



632 



635 



636 



m 637 



638 



639 



640 



0.2 



0.4 0.6 
fractional Julian day 



0.8 



Figure 1. The reduced XCOV15 time-series photometry showing the data quality and the extended coverage provided by the multi-site 
campaign. Each panel represents a whole Julian day with the actual day determined by adding 2450000 to the number on the right of 
each panel. The vertical axes are in units of millimodulation intensity [mmi] whereby 10 mmi corresponds to a 1% flux variation from 
the local mean value. The (blue) values in each panel centred around h UT represent SAAO data, the (red) values centred around 12 h 
UT represent M JUO data, while the CTIO data (light blue) is centred around 6 h UT in the 4th panel from the top and the OPD data 
(green) segments appear in panels 2,3 and 5 in the vicinity of 6 h UT. Note that the two OPD runs ra411 and ra413 have not been 
included in the plot, as they were not included in the final analysis (see discussion in the text). Also, see the online journal article for a 
colour version of this figure. 



Table 2. Journal of observations of time-series photometry of 
EC 20058-5234 obtained at Mt John Observatory in June 2004. 
Columns 3 and 6 give the UT and (modified) barycentric Julian 
date (BJD~) start times for each run and columns 4 and 5 provide 
the number of useful 10 s integrations and length of run in hours, 
respectively. 
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This is achieved by using a dichroic filter in the comparison 
star's optical train so that the red component of its spec- 
trum is transmitted to a small CCD and the remaining blue 
component is reflected to the photomultiplier tube. Contin- 
uous remote guiding or autoguiding is then possible using 
only the capabilities of the photometer. 

All of the 2004 Mt John photometry was obtained using 



the three channel photometer operated in autoguiding mode. 
A good run of clear weather in June led to a 24% duty cycle 
out of a total observing period of 7.2 days. The results of 
this work are summarised in Table [5] 



3 DATA REDUCTION 

We wish to identify periodicities in the light curve that re- 
sult from stellar pulsation. The first stage in this process is 
to convert the time series data to a form that emphasises the 
intrinsic stellar intensity changes: this requires the removal 
of a number of artifacts that result from the observing pro- 
cess. These reductions steps are relat ively straight forward 
(e.g. iNather et allll990l : lKepler|[l99ot l. but in the interests 
of completeness, and also because we have employed a novel 
smoothing procedure, we will briefly summarise our steps 
here. 

The EC 20058 field of view has two faint companions 
separated from the target by 2 and 4 arcseconds, respec- 
tively. For the aperture photometry employed in the work 
presented here, it was impractical to attempt to separate the 
target stars from its two companions (especially in indiffer- 
ent seeing conditions), so all observers used aperture sizes 
large enough to include all 3 stars and prevent any errors in- 
troduced by imperfect tracking. It is relevant to note at this 
point that the light curve pulsation amplitudes presented 
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later have not been corrected for this co mpanion star con - 
tamination, which has been estimated bv lKoen et al.l (|l995T ) 
to require multiplication by a factor ~1.4. 

For each light curve, the sky background (approaching 
50% in some cases) for each integration was estimated and 
then subtracted point by point from the raw light curves. 
For both the single channel and two channel photometers, 
sky measurements were made by moving the telescope to a 
blank section of sky at the beginning and end of each run, as 
well as several times per hour. A sky background time series 
curve was constructed for each light curve by using either 
linear or cubic spline interpolation between the measured 
sky points and then subtracting this computed curve from 
the raw light curve. 

For the three channel photometers (CTIO in the WET 
run and all of the 2004 Mt John photometry), the relative 
channel sensitivities were calibrated by measuring sky values 
for all 3 channels at the beginning and end of each run. Ap- 
propriately adjusted sky values were then subtracted point 
by point from the other two raw light curves. 

The effects of airmass changes were accounted for by 
fitting low-order polynomials to the sky-corrected data, and 
then determining fractional deviations of the data points 
from the fitted values. Although the impact of airmass 
changes can be readily modelled using a plane-parallel at- 
mosphere approximation, it is just as effective to use the 
polynomial fitting method. 

Accounting for other light-curve artifacts such as the 
effect of cloud is more problematical, but can be achieved 
for the two and three channel photometer data by normalis- 
ing the sky-corrected target data to the sky-corrected com- 
parison data for those regions affected by cloud. The effect 
of thin cirrus cloud on a moonless night up to a reduction 
amount of about 30% can normally be corrected for. This 
applies in particular to data acquired with the three chan- 
nel photometers (such as that obtained at Mt John in 2004) 
since (a) all three channels employ identical blue-sensitive 
photomultiplier tubes (Hamamatsu R647-04), and (b) the 
dichroic optical element in the comparison channel makes 
its passband effectively similar to that of the target chan- 
nel (when viewing the hot white dwarfs) , irrespective of the 
spectral type of the comparison star. However, the overrid- 
ing criterion employed here is that if a data segment appears 
irreparably cloud-contaminated then it is not included in the 
combined data set. The data sets listed in Tables [1] and [2] 
adhere strictly to this criterion. 

Any residual 'long period' variations in each reduced 
light curve segment (now expressed as fractional deviations 
from the local mean) were removed using techniques in- 
corporated in the program developed by one of us 
(DJS). This program allows the user to fit cubic splines to 
the data by marking a displayed plot with points at times 
and flux values selected by the user. The flux fitting values 
can be determined either by using a local flux average value, 
or by direct visual selection. Cubic spline fits to these points 
are then used to smooth the light curve over the entire data 
segment time interval. 

Undoubtedly, there is an increased component of sub- 
jectivity in this procedure when compared with an alterna- 
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tive more objective technique (such as high pass filtering in 
the frequency domain); the latter can be readily replicated 
by others. One might liken it to 'chi-by eye' visual fitting of 
an elementary function to a data set compared to use of a 
more objective fitting procedure, such as least squares. Still, 
the power of an experienced observer viewing data presented 
in suitable graphical form and making informed judgements 
about the quality should not be underestimated; experience 
has shown that these smoothing procedures are very useful 
in removing light-curve artifacts that are clearly not related 
to the WD pulsations of interest. 

However, the procedure is obviously open to 'abuse' 
and consequential loss (or gain) of signal. An important op- 
tion, which acts as a safeguard, allows the user to calculate 
Fourier transforms of both the input and modified data and 
view a comparison plot in order to monitor directly in the 
frequency domain the changes that are made in the time 
domain. The procedures can be easily modified or restarted 
if one is suspicious of the changes made to the light curve. 
This program has proved to be very effective in eliminating 
in the time domain obvious artifacts that are not related 
to white dwarf pulsation, such as inadequate twilight sky 
correction, and uncorrected transparency variations due to 
both airmass changes and cloud. 

The resulting corrected light curves are then, of course, 
blind to periodic variations larger than a certain period 
limit, but it is preferable to remove extraneous power in 
the time domain rather than simply rely on signal orthog- 
onality in the Fourier domain to differentiate between the 
signals of interest and other artifacts. Also, given that the 
white dwarf pulsations of interest are in the range of about 
100 s to 1000 s, then if one only makes careful use of the 
program's capabilities on longer data segments (at least an 
hour or more), then there is a negligible chance of removing 
any r eal signal. Besides, previous work (|Breger et alj|l995l . 
1996) has demonstrated that use of the continuous mon- 
itoring time-series photometry reported here is inferior to 
the more traditional 'three star' photometry when studying 
longer period (~ hour or more) variables, such as the S Scuti 
stars. 

Finally, the observation start times for each observatory 
significant data set were converted to the uniform barycen- 
tric Julian day (BJD) timescale, which corresponds to in- 
ternational time (TAI) in units of days transformed to the 
centre of mass of the solar system (jAudoin fc Guinotll200ll ; 
IStandislj[l99S| y 

ft is relevant to note here that individual integration 
times within each listed data segment were not transformed 
to the BJD timescale. This results in a maximum timing 
discrepancy of about a second (1.2 s for the 7.57 h 2004 Mt 
John June 11 data), so this omission will have a negligible 
impact on the Fourier analysis of the 10 s integration time- 
series data. However, these timescale changes are taken into 
account in software we have developed to search for long 
term period changes. 



4 FOURIER ANALYSIS 

The standard way to identify the frequency structure in a 
reduced light curve when the phase information is not re- 
quired is to calculate a power spectrum. This is presented 
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Figure 2. An amplitude periodogram (DFT) of the XCOV15 data set together with the DFT window (plot insert) using the same 
horizontal scale as the main plots. The vertical axes employ the (linear) unit of millimodulation amplitude [mma] in which 10 mma 
corresponds to a 1% amplitude modulation of the light flux. This plot makes clear that the EC20058 is a multi-periodic low amplitude 
pulsator as the two dominant modes have amplitides less than 10 mma. Readily identifiable periodicities (15 in all) are labelled and there 
is a suggestion of more real power above 6500/^Hz other than the 4 labelled peaks fi3, fi4, fn and fig. Note that £■? is readily identified 
either by comparing the structure in the DFT at frequencies just above the main f§ peak with the window function, or (better still) by 
inspecting a suitably prewhitened DFT (Fig. 3). 



for the WET data in Fig. [2] in the form of an amplitude 'pe- 
riodogram' covering the frequency range from 100 to 10,100 
nHz, and in which the vertical axes employ the scale mil- 
limodulation amplitude [mma]. We have employed an am- 
plitude scale rather than a power scale (which is now com- 
mon practice) to represent signal 'power', as this facilitates a 
more direct intuitive connection with actual sinusoidal vari- 
ations in the light curve and comparisons with least squares 
fitting of sinusoidal functions. Following the work of lScargQ 
(1982), a number of authors use the term 'amplitude peri- 
odogram' or simply 'periodogram' for these plots. We will 
henceforth simply use the acronym DFT (for discrete Fourier 
transform) to indicate this procedure. 

The WET data cover a period of 8.3 days, which cor- 
responds to a frequency resolution of ~ 1.4 fiHz. The insert 
panel in Fig. [2] illustrates this resolution by plotting the 
DFT 'window' for the transform using the same horizontal 
scale width as the main plot. This window corresponds to 
the DFT of a synthetic noise-free sinusoid calculated at the 
same times as all the light curve data: it directly illustrates 
the frequency resolution and the impact of the spurious side 
lobes, or 'aliases', that result from turning the signal on and 
off during the observing period. The central plot value has 



been shifted to zero frequency for clarity. The near contin- 
uous WET data has minimal alias interference, which is, 
of course, the primary reason for the observation strategy 
pursued by the WET collaboration. Note that in the wider 
Fourier literature, the additional peaks in a DFT deriving 
from a particular observing window are termed 'spectral 
leakage', while the term 'alias' is reserved for effects near 
the Nyquist frequency resulting from undersampled data. 
Since we are considering frequencies in our well-sampled 
data much smaller than the Nyquist frequency (50,000 /dTz), 
following other authors, we will use the term 'alias' for this 
phenomenon, as it more accurately describes the additional 
forest of sometimes confusing peaks that occur, particularly 
for extended data sets with regular daily gaps (eg, see the 
right hand panels in Figs. [5] and [6] corresponding to the 
MJ04 data). 

Direct visual inspection of the DFT in Fig. [2] allows 
one to readily conclude that there are about fifteen 'real' 
periodicities present in the data. To this end, the fifteen 
visually obvious frequencies, ranging from ~ 2000/iHz to 
10 4 /dTz, have been labelled fi to fis in order in the plot, and 
also listed in Table [3] At this stage a somewhat conservative 
(but arbitrary) threshhold has been used to decide that the 
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Figure 3. The DFT of the XCOV15 light curve after prewhitening (see text) by the 10 periodicities with the highest amplitudes: these 
frequencies are indicated by the downward (red) arrows. The vertical axes use the same units as in plot 2 and the horizontal dashed 
lines at 0.54 mma in the plots correspond to the 0.001 false alarm probability detection threshhold established by the Monte Carlo 
data-shuffling method discussed in the text. Using this value we can assert that there is real power at ~ 4800 /iHz (fg), but not at the 
frequency just below 4000 /iHz marked a. The vertical (green) lines in the lower panel, annotated with the letters b - h, correspond to 
the predicted positions of combination frequencies where there are indications of real power in the DFT (see Table [3J. The detected 
frequencies labelled fi5 and fi6 having amplitudes just above the significance threshhold also correspond to the predicted values of 
combination frequencies. 



peaks fi3 and fis, for example, are real and other smaller 
amplitude peaks (such as fg and others at higher frequencies) 
are uncertain. Note that the vertical amplitude scales in the 
upper and lower panels of Fig. [5] are different, and that the 
reality of {? is clear from either a careful comparison of the 
DFT in the region of f$ and the form of the window function, 
or from the analysis presented below. 

The EC 20058 light curve is dominated by the beating 
between the two largest ~ 8 mma pulsation modes at 281 s 
and 257s (f6 and is, respectively). With these modes com- 
bining to produce a maximum modulation of less than 2% (~ 
17 mma) every 50 minutes, this star is certainly not a large 
amplitude pulsator. However, as noted previously, the ac- 
tual intrinsic stellar flux amplitudes are approximately 40% 
larger due to the uncorrected companion star contamina- 
tion. 



4.1 Light curve prewhitening 

A very useful procedure for separating peaks in the DFT 
corresponding to real power from those peaks produced by 
window effects and/or noise is that of 'prewhitening'. In this 



procedure, one or more selected frequencies are first removed 
from the light curve by subtracting least-squares fitted si- 
nusoids, and then a DFT of the 'prewhitened' light curve 
is computed. This method is particularly enlightening when 
one is attempting to identify small amplitude periods that 
are mixed in with the alias pea ks of a much la rger amplitude 
perodicity. As pointed out by IScarglel l|l982h . least-squares 
fitting of sinusoids in the time domain and computing DFTs 
are essentially identical analytical procedures, so they both 
should independently lead to the same conclusions. Conse- 
quently, a careful comparison of the window function with 
the DFT spectrum in the vicinity of a peak should provide 
the same information as that revealed by the prewhitening 
procedure. However, uncertainties in the plots sometimes 
dominate such comparisons, and therefore the prewhitening 
method is often more decisive, as some examples below make 
clear. 

Fig. [3] shows the WET DFT covering the same fre- 
quency range as for Fig. [2] but with the light curve first 
prewhitened by the 10 highest amplitude frequencies: fi - 
f4, f6, fg, fio - fi2 and fi4. The y axis scales in this plot have 
been expanded, and the removed periodicities are indicated 
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Table 3. Detected frequencies in the DFTs for three EC 20058 
data sets. The amplitudes (in millimodulation units) for the 1997 
XCOV15 WET run are given in column 5, while the corresponding 
values for the 2004 Mt John data and the SAAO 1994 discovery 
data are given in columns 6 and 4, respectively. The adopted com- 
bination frequencies are listed in the last column and amplitude 
values below the relevant significant levels are given in parenthe- 
ses. Note that the amplitude values have not been adjusted to 
account for the companion field star light contamination and the 
1994 values have been obtained from a subset of the Koen et al. 
(1995) data - 4 runs in July 1994. 



item 


freq. 


period 


amplitude [ 


mma] 


combn. 




[/iHz] 


N 


1994 


1997 


2004 


freqs. 


fi 


1852.6 


539.8 




2.0 






h 


1903.5 


525.4 




1.9 


1.5 




h 


2852.4 


350.6 


(1.6) 


1.4 


1.3 




h 


2998.7 


333.5 


2.2 


3.0 


2.6 




f s 


3489.0 


286.6 




1.7 


1.3 




fe 


3559.0 


281.0 


8.8 


8.5 


7.6 




h 


3640.1 


274.7 




1.0 






h 


3893.2 


256.9 


8.0 


8.4 


7.4 




a 


3924.2 


254.8 




(0.5) 






fg 


4816.8 


207.6 




0.7 






fio 


4887.8 


204.6 


1.2 


2.7 


2.4 




t'n 


4902.2 


204.0 




1.7 


1.4 


f 2 +f 4 


fl2 


5128.6 


195.0 


3.5 


2.5 


1.8 




b 


5462.5 


183.1 




(0.5) 




h+h 


c 


5745.8 


174.0 




(0.3) 




fi+fs 


fl3 


6557.6 


152.5 




0.6 




f 4 +f 6 


d 


6745.6 


148.2 




(0.5) 




f 3 +fs 


c 


7118.0 


140.5 




(0.5) 




2f 6 


hi 


7452.2 


134.2 


2.3 


1.7 


2.0 


fe+fs 


f 


7533.5 


132.7 




(0.4) 




f 7 +fs 


g 


8127.3 
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(0.5) 




f 4 +fl2 


h 


8446.8 


118.4 




(0.5) 




fe+fio 


fir, 


8461.2 


118.2 




0.6 




f2+f 4 +fe 


fl6 


8781.0 


113.9 




0.6 




fs+fio 


fl7 


9021.7 


110.8 


1.6 


0.9 


(1.1) 


f 8 +fl2 


fl8 


10030.7 


99.7 




0.7 




f 2 +f 4 +fl2 



by the downward arrows. The prewhitening was carried out 
by simultaneously fitting (via least squares) 10 sinusoids to 
the data using periods obtained from the original DFT. If 
the periods are kept fixed, linear least squares fitting can be 
used, and in practice this is adequate for the well-defined fre- 
quencies in the data (especially if accompanied by a limited 
grid search). 

An examination of the Fig. [3] DFT shows that there 
are two clearly significant peaks (fs and ij) on either side of 
the main peak (fe). In addition, there are two peaks in the 
upper panel (a and fg) that look as though they could be 
significant, a collection of 9 peaks in the lower panel (labelled 
b - h, fis, fig), many of of whom could be significant, and one 
or two peaks at low frequencies that need to be investigated. 
These latter peaks will be discussed in the next section, but 
clearly we need a quantitative criterion for distinguishing 
between real and noise peaks. 

A detailed inspection of the prewhitened DFT for the 
entire data set listed in Table Q] in the vicinity of the dom- 
inant periods % and fs revealed that there was a small but 
significant amount of power left in the form of residual 'win- 



dow mounds' for both frequencies. In other words, the sub- 
traction of the two periodicities from the light curve had not 
removed all of the power at those frequencies. For a pulsat- 
ing white dwarf exhibiting very stable and highly coherent 
luminosity variations, this can be caused by one or more 
data segments having a timing error. By a process of trial 
and error, the 'culprits' were determined to be the two OPD 
observation files ra411 and ra413 marked with an asterisk 
in Table Q] Upon deletion of these two files from the over- 
all light curve, the prewhitened DFT produced the results 
graphed in Fig. [3] and shown in more detail in Fig. [5] (left 
panels): no signatures of the subtracted frequencies remain. 
This discovery illustrates the utility of the prewhitening pro- 
cedure, as the 'contribution' from the two relatively small 
data segments could easily have been overlooked otherwise. 
In order to be certain, the timing discrepancies were inves- 
tigated graphically by overlaying a time series plot of the 
two data segments with the predicted light curve obtained 
from a least squares fit of the two dominant frequencies to 
the rest of the WET data. Both plots clearly revealed the 
presence of a timing error. 

Given the unlikely nature of the alternative hypothesis 
— the star was behaving differently during these periods — 
subsequent analysis removed these data from the combined 
WET light curve. No obvious explanation for the timing 
error has been found. 



4.2 DFT noise simulations 

A useful quantitative criterion for differentiating between 
signal and noise peaks in the DFT is the concept of th e 
false alarm probability (FAP) introduced bv lScargiel l|l982t ). 
which makes the statistical nature of the process explicit. 
Using this concept, one computes the probability of a given 
peak in the DFT being due to noise for some chosen thresh- 
hold. A high threshhold for false positives can be set by 
ensuring that this probability is small, and only counting 
peaks that are above this threshhold as real. 

Although the FAP can be estimated theoretically on the 
basis of some ass umed model for the noise characteristics 
(e.g. IScarglelll982l ). it is preferable and more representative 
of the actual data to determine it by a Monte Carlo simu- 
lation method. And, assisted by the speed of modern com- 
puters, this is a practical proposition for even large data 
sets. Such a simulation has been carried out for the data 
sets discussed here using the program tsmrarjf] developed by 
one of us (DJS). The results of this work for the WET data 
set are that peaks with an amplitude of 0.54 mma have a 
FAP value of 0.001, meaning that there is only one chance 
in 1000 that such a peak can be produced by a random noise 
'conspiracy'. The corresponding value for the Mt John 2004 
(MJ04) data set is 1.2 mma. 

These threshholds are marked as horizontal dashed 
lines in appropriate parts of Figs. [3] [5] and [6] showing the 
prewhitened DFTs for both data sets. 

Although the simulation procedure is conceptually 
straight forward, there are several points that should be 
noted and a number of organisational tools are required for 
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Figure 4. Histograms showing the results of the Monte Carlo data shuffling exercises (see text) invoked to establish a false alarm 
probability detection threshhold (0.001) for real power in the prewhitened DFTs. The horizontal axis gives the magnitude of the 
maximum peak in each DFT of the randomised time series data jand the vertical axis gives the number of occurrences per bin interval. 
The two histograms (blue and green) centred around 0.4 mma correspond to separate 1000 shuffling trials of the pre-whitened XCOV15 
data and use a bin interval of 0.005 mma, while the histogram on the right employs a bin interval of 0.01 and represents the same process 
applied to the 2004 multi-night Mt John data. For each DFT of the randomised time series, the maximum peak A m in the range 1600 
- 10,000 /iHz was selected. The downward arrows mark the maximum peaks obtained in all 1000 trials (A™ ax ) and 4 times the average 
DFT peak height (see text) for each ensemble of the DFTs (4 <A av >)■ Sec the online journal article for a colour version of this figure. 



a flexible implementation. So, we briefly outline the tsmran 
methodology here. 

The essential idea is to start with the reduced light 
curve that has been prewhitened by the 'clearly real' fre- 
quencies, randomly rearrange the time order of the data 
points, compute a DFT of these randomised data covering 
the frequency range of interest, and then record both the 
maximum (A m ) and the average (A av ) value of peak heights 
in the DFT. Randomising the time order of the data destroys 
the coherency of any periodic signal remaining in the light 
curve but preserves the uncorrelated noise characteristics of 
the data. The highest peak in the DFT of a time-shuffled 
data set is then a one sample estimate of the maximum ex- 
cursion that can occur due simply to random noise effects. 
If this procedure is repeated a large number of times (say 
1000), and the highest peak in each DFT recorded, then one 
can infer that the maximum value, A™ ax , in this ensemble 
of highest peak values is a direct estimate of the DFT am- 
plitude threshhold for a false alarm probability at the 1 in 
1000 level. 

It is instructive to plot a histogram of the recorded max- 
imum peaks (A m ) for each DFT: 1000 member ensembles 
for both the WET and MJ04 data are depicted in Fig. g| 
There are two separate such histograms (blue and green) 
for the WET data centred around 0.4 mma, and one for the 
MJ04 data centred around 0.9 mma. It is clear that the less 
comprehensive MJ04 data exhibits larger A m values, as one 
would expect. 

The program tsmran also records the average peak 
height for each randomised DFT, partly as a check on its 
operation, and as one might expect there is little variation 
from DFT to DFT. The histograms of these average peak 
height values have not been shown, but four times the en- 
semble average (4 <A av >) for each data set has been indi- 
cated in Fig. [4] as one sometimes sees a criterion similar to 



this u sed in the literature to identify real power (e.g. iKeplerl 
1993). Note, however, that we are using here (the average 
of) the average amplitude DFT values and not the square 
root of the average power spectrum values. 

A few more comments about the algorithms employed 
in tsmran are relevant. For the relatively large data sets 
considered here (WET ~ 45k and MJ04 ~ 15k integrations, 
respectively) a long period random number generator is re- 
quired: the routin e ran2 from the Numerical Recipes suite 
Fress et alJll992h was found to be equal to the task. One 
also makes significant efficiency gains by using a fast Fourier 
transform (FFT) algorithm to compute the DFTs f or suit- 
able data sets; the routine realft from IPress et al.l accom- 
plishes this task in tsmran. 

However, since FFT algorithms assume that the time- 
series data correspond to equal contiguous sample intervals, 
some preliminary organisation of the data segments is re- 
quired before a FFT can be employed. Since all the ob- 
servations reported in this paper correspond to 10 s inte- 
grations it was possible to create an array of integrations 
that satisfied this attribute (with adequate precision for the 
intended purpose) without resorting to such (undesirable) 
procedures as interpolation. All deleted integrations in each 
'significant' data segment were restored using zero padding, 
and the shape of the overall DFT window function was es- 
sentially preserved by separately randomising the points in 
each segment before combining them together in one overall 
time series, all the while maintaining the data gaps using 
zero padding. Last, the combined time series is extended to 
an integer power of two (as required by realft) using zero 
padding, in order to maximise the efficiency of the FFT 
procedure. 

A comment on the data timescale is appropriate. Al- 
though the above assumption of uniform contiguous inte- 
gration intervals represents the original terrestrial TAI ob- 
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Figure 5. Expanded amplitude DFTs of the frequency region near modes fg and fg, comparing the 1997 XCOV15 data set (top left) with 
the June 2004 (Mt John) data set (top right). The panels at the bottom, using an expanded amplitude scale, depict the corresponding 
DFTs of each time series prewhitened by the frequencies fg and fg (blue curves). The periodicities corresponding to fs and {j are clearly 
evident in the prewhitened XCOV15 data (bottom left) but only fs is evident in the Mt John data. There is a suggestion of power (marked 
'a') in the XCOV15 data adjacent to fg that is just below the adopted detection threshhold (dashed line). Both bottom panels also include 
overlay plots (in red) of prewhitened DFTs in which the remaining periodicties have been removed. See the online journal article for a 
colour version of this figure. 



servation timecale, offsets are introduced after the transfor- 
mations to the BJD timescale. This will have little impact 
on use of the FFT for the simulation exercise, and negligi- 
ble impact anyhow if the offsets are not very large. For the 
record, the total offset change across the WET data is ~ 7 s, 
and for the MJ04 data it is ~ 25 s. 

Individual data segments corresponding to each run are 
read separately into tsmran using a list of file names in an 
input file provided by the user. A glance at Fig. [T] (especially 
the colour version) shows that there is some data segment 
overlap (eg, the fourth panel from the top: SAAO, CTIO and 
MJUO on 6 July). These overlaps were removed by deleting 
the minimum number of points in order to allow contiguous 
data to be created. 

A version of the program exists which uses the much 
slower (computationally inefficient) direct DFT algorithm. 
This then obviates the need for much of the data organi- 
sation code in tsmran, which admittedly was a significant 
proportion of the programming coding effort. This method 
is more general and can be invoked for non-uniform time 
series data sets, but it is significantly slower. In contrast, 
the FFT method is extremely fast: on a 'standard' modern 



laptop the one thousand 64k WET DFTs are computed in 
less than a minute — effectively 'on the fly', one might say. 

At this point it is relevant to emphasise that our sim- 
ulations are dealing with a random noise model only; other 
non-random sources of 'noise' (eg, periodic telescope drive 
error) need to be handled on a case-by-case basis. In addi- 
tion, we limited our region of interest to frequencies above 
1600^iHz, as (a) there are no obvious peaks in the DFT in 
this region that stand out above the forest of other peaks and 
(b) we would need to consider a higher detection threshhold 
due to the increasing impact of residual sky noise. 



4.3 Detected frequencies 

Referring to Figs. [21 El [S] and |S] and Table [3] we have iden- 
tified 18 definite frequencies (f*) in the WET data as being 
real, using the adopted 0.001 false alarm probability. The 
amplitudes for these frequencies are given in column 5 of 
the table, while the amplitudes of other possible frequency 
detections (a - h) are listed in the table within parentheses. 

Two frequencies that merit further discussion are 'a' 
and and fg; they are depicted in more detail in the panels of 
Figs. [5] and [6] and . These frequencies are essentially outside 
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Figure 6. Expanded amplitude DFTs of the frequency region near modes fio and fu, comparing the 1997 XCOV15 data set (top left) 
with the June 2004 (Mt John) data set (top right). The panels at the bottom depict the corresponding DFTs (blue curves) of each time 
series prewhitened by the frequencies fio and fu. The horizontal dashed lines in the bottom panels represent the detection threshholds 
for significant power in the respective data sets (see discussion in Section 4.3). The fg frequency is above the detection threshhold in the 
WET data and, although it is below the Mt John data detection threshhold, it does appear to make its presence felt. The red overlay 
plots in both bottom panels represent DFTs of each data set further prewhitened by the f9 periodicity. See the online journal article for 
a colour version of this figure. 



the region where the issue of linear combinations arises (see 
next section) and their amplitudes are very close to the 0.54 
mma significance threshhold. They only first appear as in- 
teresting in the prehwhitened DFT of Fig El and the MJ04 
data set provides additional information. This is depicted in 
the right panels of Figs. [S]and[5]by graphing both the MJ04 
DFTs and prewhitened DFTs covering the same frequency 
ranges as for the WET DFTs. 

The XCOV15 case for fg is clearly presented in the two 
left panels of Fig. [6] The top left panel shows the DFT with 
this frequency marked (along with fio and fu), the lower 
left panel (blue curve) shows the DFT prewhitened by fio 
and fu and power at the fg frequency above the 0.54 mma 
threshhold, and the dark (red) curve shows the DFT fur- 
ther prewhitened by the fg frequency in which the peak has 
disappeared. 

The two right panels in Fig. [6] present the MJ04 data 
case in the same format as the left panels. There is clearly 
evidence of power in these data at the frequency correspond- 
ing to fg which follows the shape of the window function and 
whose maximum peak height is not far below the relevant 1.2 
mma significance threshhold. Also, the further prewhitening 
procedure removes these peaks from the DFT (red curve). 



Although the principal peak is below the 1.2 mma MJ04 
data significance threshhold, so in isolation would not be in- 
terpreted as a detection using our stated threshhold, it does 
provide definite supporting evidence for the XCOV15 detec- 
tion. 

On the other hand, the case for (and against) 'a' is pre- 
sented in a similar manner in Fig. [5] In the lower left panel, 
the blue curve depicts the XCOV15 data prewhitened by the 
two dominant frequencies (fe and fg) and the expanded am- 
plitude scale clearly shows the reality of fs and {7, and also 
a possibly real frequency 'a', tantalizingly just below the 
0.54 mma significance level. The dark (red) curve is a DFT 
further prewhitened by these three frequencies and empha- 
sises the point that the 'a' peak could represent real power. 
However, the same procedure for the MJ04 data tells a dif- 
ferent story. The blue curve in the lower right panel shows 
the MJ04 DFT prewhitened by the dominant f6 and fg pe- 
riodicities and there is no evidence of 'a'. 

It is interesting that there is also no evidence for i~j in 
the MJ04 data. However, there is without question power at 
this frequency in the WET data so further confirmation is 
not required. It does illustrate, though, that power in this 
quite stable pulsator at various frequencies is still variable at 
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a low level, the result of real physical amplitude instabilities, 
or perhaps beating. 

In view of the above, we have included f9 as a real de- 
tected periodicity, but not included 'a' — we leave it in the 
summary table as a bracketed 'maybe'. 



4.4 Linear combination frequencies 

Since the mechanism that converts mechanical movement of 
stellar material to luminosity variations at the surface ap- 
pears to be nonlinear, both harmonics and combinations of 
the basic mode frequencies can be expected to appear in the 
light curve. At least, this is the conclusion of many previous 
studies of the pulsating white dwarfs — especi ally the large 
ampl itude pulsators, such as DBV GD 358 |Kepler et al.l 
12003! ). 

The last column in Table [3] lists whether any of the de- 
tected frequencies in the light curve can be identified as a 
linear combination of other detected frequencies, as well as 
identifying whether possibe linear combinations could ac- 
count for 'suspiciously' large peaks that are below the 0.54 
mma detection threshhold; all these quantities are also in- 
dicated on the plot in Fig. [3] 

So, all of the actual 7 detected frequencies fn and fi3 - 
fis can be explained as combination frequencies. This inter- 
pretation leaves fi2 as the highest frequency pulsation mode. 

The two largest combination frequencies with ampli- 
tudes ~ 1.7 mma merit further discussion. The frequency 
fi4, being simply the sum frequency of the two highest am- 
plitude modes (f6 and fg) is obviously the first combination 
to look for, but fn being considered a combination of the 
lower amplitude modes f2 (1.9 mma) and f4 (3.0 mma) is 
perhaps surprising. Furthermore, these two modes appear to 
combine separately with f6 and fi2 to produce signal power 
at fis and fig, respectively. 

Whereas some of the higher frequency combinations 
would be difficult to explain as direct mode frequencies us- 
ing realistic white dwarf models, this is not true for the fre- 
quency fn: it is not absolutely certain that this frequency 
results from nonlinear combination effects. However, given 
the frequency matches mentioned above, and the likely dif- 
ficulty of explaining the presence of the four closely spaced 
frequencies fg - fi2 in terms of low-order £ pulsation modes 
(see Fig. [7] and discussion next section), we adopt the con- 
servative approach and delete fn from the list of inferred 
modes. 

Column 5 in Table lists the 11 frequencies (fi - fio 
and fi2) that we consider represent pulsation modes from 
the work presented here, while columns 4 and 6 in the same 
Table li st for compa r ison a ll the frequencies detected in the 
work of iKoen et all (|l995l ) and independently in the MJ04 
data set, respectively. 



are, in fact, two of these: one is a sequence of periods with 
approximately equal spacing (which is obviously best viewed 
in period space), while the other is 'mode-splitting', which is 
seen more clearly in frequency space. Second, we can adopt 
a global fitting procedure which endeavours to match the 
predicted modes to the data, su ch as that pioneered by 
iMetcalfe. Nather fc Wingetl (|2000l 'l. 

We need two basic results from the relevant pulsa- 
tion theory in order to look for any systematic trends. The 
buoyancy-driven g-mode pulsations in white dwarfs (e.g. 
lUnno et all 1 19891 ; Hansen et all [20041 ) can be characterised 
by three integer quantum indices: n (or sometimes k), £ and 
m. These indices specify a unique eigenmode of oscillation 
in which n determines the radial form of the eigenfunction 
and £ and m specify the angular behaviour via a spherical 
harmonic function Y^ in (9, <f). If spherical symmetry holds 
(no rotation and/or magnetic field), then the mode periods, 
Yl n i, are independent of the index m and in the limit of large 
n we can write 



n,- ( 



IT; 



n + e 



< AIT > : 



IT, 



The quantity IIo is inversely proportional to an integral of 
the Brunt- Vaisala (or buoyancy) frequency over the star pro- 
file, and hence decreases with increasing stellar mass (due 
to increasing local gravitational field strengths). Hence, Ho 
establishes a basic timescale for the period structure; e can 
be neglected for our purposes. 

Rotation of the star removes the m-degeneracy of the 
mode frequencies, and to first order leads to mode frequen- 
cies that are a function of all three quantum indices in the 
form 



V n tm — Vnt + W,(l — C n e)Cl 

where u n i is the degenerate frequency, m takes the 2^ + 1 in- 
teger values between — £ and +£, Q is the rotation frequency 
of the star and C n e takes the approximate dimensionless 
form C n i = l/£(£ + 1) for high radial overtone modes. The 
net effect of this is that all (n,£) g-mode pulsations undergo a 
'Zeeman-like' splitting due to the rotation of the star which 
can theoretica lly be seen as 2£+l closely spaced frequencies 
fsee [Cro3[l98i for a nice physical explanation). 

A magnetic field also destroys the spherical symmetry 
and leads to frequency splitting: in the case of a small mag- 
netic fi eld aligned with the pulsation axis, £ + 1 splitting 
occurs jjones et al.lll989h . Whether any of these modes are 
sufficiently excited such that they are detectable is another 
story. 



5 COMPARISON WITH MODELS 

Our aim here is to compare the 11 detected pulsation modes 
for EC 20058 with predictions from white dwarf models in 
order to tell us something about the star. There are two 
basic approaches to this task. In the first instance we can 
look for relatively simple systematic trends in the observed 
pulsation spectrum that depend on stellar properties. There 



5.1 Model calculations 

In order to provide concrete examples, we have determined 
a range of pulsation modes for a sequence of WD models 
with effective temperatures between 30 000 and 20 000 K. 
The results of our calculations are displayed in Fig. [7] (left 
panel), along with the detected mode periods for EC 20058 
(middle panel) and the DBV class protype (GD 358, right 
panel) . 
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Figure 7. A pulsation mode diagnostic diagram showing the predicted periods of the pulsation modes as a function of effective temper- 
ature for a range of 0.6Mq white dwarf models. The left panel graphs all I = 1 modes (red solid lines) in the period range and only some 
of the 1 = 2 modes (blue dash-dot lines). The middle panel shows the detected modes for EC 20058 and the right panel shows the modes 
detected for GD358 for comparison purposes. See text for more details and the online journal article for a colour version of this figure. 



The WD models and the respective pulsation frequen- 
cies were calculated using code that has itself evolved over 
the years, and is described in a succe ssion of PhD theses at 
the University of Texas, Austin - see iMetcalfe et al.l (|2000l ) 
for a brief description and appropriate references. For our 
fiducial model, we chose a mass of 0.6 Mq, a 50/50 C/O 
core and a fractional helium layer mass of 10~ 3 . Six evolved 
WD models with effective temperatures at 2000 K inter- 
vals between 30 000 K and 20 000 K were computed, and the 
pulsation modes determined for each of these. The pulsa- 
tion modes for models of intermediate temperatures between 
these values were estimated by fitting cubic splines to all the 
mode periods as a function of effective temperature and the 
results are plotted in Fig. [7] all the I = 1 mode periods 



are represented by solid (red) lines in the displayed period 
range, and a few of the shorter period 1 = 2 modes are dis- 
played using the (blue) dash-dot lines. Since we are using 
these model calculations to aid the phenomenology discus- 
sion below, the precise values of the model parameters are 
unimportant. 

A very obvious feature of the pulsation period 
structure in the WD models is that all mode peri- 
ods increase with decreasing effective temperature (e.g. 
iBradlev, Winget fc Woodl I1994T I. and therefore increasing 
age. Over human timescales this effect is very small, but 
nevertheless the extremely high stability of many white 
dwarf pulsators means that this effect can be observed 
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(e.g. iKepler et all lOOSbh . We will discuss the prospects for 
EC 20058 in the last section. 

We have chosen to discuss only the lowest order I = 1 
and 1 = 2 pulsation mod es, as we make the usual assump- 
tion (lDziembowski|[l977l) that the effect of geometrical flux 
cancellation across the observable stellar disk should ensure 
that modes with larger £ values, if excited, contribute little 
to the observable flux changes. However, at this point it is 
appropriate to sound a small note of caution by mentioning 
that there is evidence of a pulsation mode in the DAV pul- 
sator P Y V ul (G 185-32) that has an £ value of at least 3 and 
possibly 4 (|Thompson et al" ] 12004 lYeates et alj|2005h . It is 
worth noting, though, that this object exhibits an unusual 
pulsation spectrum. 

For the discussion in the next section we will focus 
on some details of the 28 000 K model, since this is clos- 
est to the effective tem perature obt ained for EC 20058 b y 
iBeauchamp et all (jl999h , but see also ISullivan et aD (|2007h . 
The pulsation modes in Fig. [7] for this model exhibit a mean 
period spacing of 37.4 s and a range of 33.4 - 43.0 s for the 
I = 1 modes, and a mean value of 20.4 s with a range of 16.6 
- 23.5 s for the 1 = 1 modes. 

Among other things, these model values demonstrate 
that the pulsation theory summarised above does only pre- 
dict approximately constant period spacings for a sequence 
of modes with a given I and varying n. Also, this theory 
predicts that the ratio of the 'constant' period spacings for 
sequences of £ = 1 a nd I = 2 in the star should equal 
•y/2(2 + l)/- v /l(l + 1) = 1.73. However, this ratio for the 
mean period spacings in our computed model (covering the 
period range in Fig. is 37.4/20.4 = 1.83. This is close to 
the approximate asymptotic theory value, but not identical, 
as we are in the low n regime. 

5.2 Period and frequency phenomenology 

It follows from the previous discussion that in the first in- 
stance we should look for a sequence of observed periods 
with approximately equal spacing in the pulsation spectrum 
of EC 20058. As is clear from Fig. [7] such a pattern was de - 
tected in CD 358 dWinget et all 1 1994 IKepler et af]|2003h . 
and they were all readily interpreted as having spherical 
degree 1=1 due to the fact they displayed clear triplet (ro- 
tational) splitting. This conclusion was possible, inspite of 
the complicated nature of the CD 358 pulsation behaviour. 

No such simple pattern is immediately apparent in the 
EC 20058 spectrum. Nevertheless, we will attempt to iden- 
tify any trends. If one assumes that a number of modes in 
some sequence are not excited above an observable thresh- 
hold, then one could consider some of the pairs fi-f2 (AP = 
14.4s), f 3 -f 4 (17.1s), f 6 -fg (24.1s) [or f 7 -f 8 (17.8 s)] and f 9 - 
fi2 (12.6 s) as visible members of such a sequence. But even 
ignoring the close spacing between the pairs, the inter-pair 
spacings do not match any assumed reasonable sequence. 
Only the two pairs £3- £4 and fV-fs show a similar period 
spacing with a mean value of about 17.5 s. However, the pe- 
riod spacing between these pairs is not even close to being 
an integral multiple of this mean value. Even if it was, we 
would then face the task of interpreting the small period 
spacing in terms of either a sequence of only 1=2 excited 
modes and/or an improbably large model mass ~ O.&Mq 
(as period spacing decreases with model mass). 
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Figure 8. A schematic amplitude spectrum illustrating the pro- 
posed rotational and magnetic mode splitting model. 



Another possibility is to temporarily put aside the ques- 
tion of the model mass and consider the two pairs fe— £s 
(24.1s) and fi-f2 (14.4 s) as members of separate 1=1 
and 1 = 2 sequences, respectively. The period spacing ra- 
tio is then 24.1/14.4 = 1.67, and this is consistent with both 
the predicted pulsation theory value (1.73) and the actual 
model estimates given in the previous section. However, in 
addition to concerns about the implied total model mass, 
it is hard to argue convincingly that only two pairs of fre- 
quencies are clear evidence of the predicted sequence. Also, 
we have not independently established any I values for the 
modes, using for example rotational splitting, as was suc- 
cessfully exploited for CD 358. 

The most likely scenario is that we are seeing a com- 
bination of I = 1 and 1 = 2 excited modes, coupled with 
possible rotational frequency splitting. An inevitable con- 
clusion is that the two modes with the largest amplitude 
(fe ~ 281s and fs ~ 257 s) correspond to different £ values, 
presumably 1=1 and £ = 2 in some order. 

We now investigate the possibility of rotational splitting 
of pulsation modes in the frequency spectrum. A qualitative 
inspection of Figs. [5] and [3] suggests that some of the modes 
f&i hi H an< i even the 'not accepted' frequency 'a' might be 
caused by rotational splitting. Using the A/ (/xHz) values 
fe — fs = 70 fiHz, ft — fe = 81 fiHz, a — fs = 31 piHz and 
fio — /q = 71 fiHz, we find that the first and fourth of these 
pairs lead to similar splitting values of ~ 70 /iHz. If these 
phenomena are due to splitting, and the modes in question 
(fe ~ 281s and fio ~ 205 s) have 1 = 1, then the white 
dwarf is rotating with a period of close to 2 hours. We are 
of course assuming that only one non zero m component is 
excited to an observable amplitude out of a possible two for 
each dipole mode. Note that this interpretation also provides 
further support (if any is needed) for our argument that 
real power is present at the frequency fg. A less likely £ = 2 
interpretation for both modes would requires only one out of 
4 possible non zero m value modes in each case to be excited, 
and would predict a rotation period close to 4 hours. 

Given the deductions of other white dwarf rotation 
rates, a ~ 2 hour period represents a rapidly rotating white 
dwarf. This rotation rate is certainly not unphysical, so we 
tentatively put it forward as a possible interpretation, al- 
though the evidence is hardly overwhelming. Note that the 
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Table 4. Optimal model parameters for EC 20058. 



Parameter 


C 


C/O 


O 


Average 




(pure) 


50:50 


(pure) 




T eff (K) 


28100 


28 200 


28 400 


28 250 ±150 


AL (M Q ).... 


0.550 


0.550 


0.550 


0.550 ±0.005 


bg(Afenv/M*) 


-3.56 


-3.60 


-3.62 


-3.59 ±0.03 


log(M Hc /M*) 


-6.42 


-6.46 


-6.48 


-6.45 ±0.03 


up (sec) 


1.89 


1.94 


2.03 





pre-white dw arf PG 2131±066 ha s a measured rotation rate 
of ~ 5 hours (|Kawaler et al.lfl995T ). Interestingly, the issue of 
measured white dwarf rotation rates is somewhat controver- 
sial as large angular momentum losses via some mechanism 
are required in order to make these measured rates con- 
sistent with the known rotation rates of their much larger 
progenitors. 

Before proceeding to the global modelling analysis pre- 
sented in the next section, we will include the small (1 mma) 
satellite frequency fy in the mode splitting model that we 
adopt. A glance at Fig.0 or the schematic version in Fig. [HJ 
reveals that fs and h appear to form an asymmetric pair 
either side of the large amplitude 'parent' mode % (~ 281 s) 
with separations of 70s and 81.3s, respectively. There are 
two ways to produce asymmetric splitting: second-order ro- 
tation effects and the inclusion of a mag netic field. 

Second order rotation effects (e.g. IChlebowskil Il978l ) 
have been detecte d in the pulsation spectrum of the DAV 
pulsa tor L19-2 jO'Donoghue fc Warner] Il982l ; ISullivanl 
1 19981 ). The dominant 192 s mode for this pulsator exhibits 
two low amplitude satellite modes that are separated by ~ 
13 ^Hz from the main peak, and a rotational splitting ex- 
planation yields a rotation period of ~ 10 hours, assuming 
an I — 1 value for the mode. Furthermore, the frequency 
splitting shows a 1% asymmetry, which is consistent with 
a second order rotational effect. Even though EC 20058 is 
estimated to rotate five times faster and the second order 
effect depends on the square of the rotation frequency, this 
is not enough to explain the asymmetry displayed in Fig. [8] 

The combined effect of rotation and a magnetic field will 
also produce asymmetrical splitting. Thus, if the frequency 
separations adjacent to the fg (281 s) mode are divided into 
a ~ 75 /iHz symmetrical rotational splitting and a magnetic 
field frequency shift (increase) of ~ 5 ^iHz, then this would 
explain the observed triplet. A magnetic field of ~ 3 kG and 
the 2 h rotation period would accomplish this (IWinget et al.l 
1 1 99-1 appendix). As suggested in Fig. [8] evidence of a similar 
triplet structure around the fio mode would add consider- 
able weight to this proposal, but as there is no evidence 
of power at fio+ ~ 81 /iHz, we simply adopt these implied 
constraints in preparation for the global modelling presented 
below. 

5.3 Model period fitting 

Starting with a list of detected pulsation modes and then 
endeavouring to deduce the parameters of a specific stellar 
model that reproduces these modes is a classic example of 
an inversion problem in physics. These tasks are notoriously 
difficult and are plagued by the lack of unique solutions. 
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Table 5. Periods and Mode Identification. 



^calc 

fnum Pobs Pure C 50:50 C/O Pure O t n 



fl2 


195.0 


192.57 


192.82 


192.76 


2 


6 


fio 


204.6 


205.23 


205.72 


205.86 


1 


3 


f 8 


256.9 


260.37 


260.61 


260.85 


2 


9 


fe 


281.0 


280.84 


281.08 


281.13 


1 


5 


U 


333.5 


330.89 


330.65 


330.70 


2 


12 


h 


350.6 


350.44 


350.44 


350.39 


2.1 


13,7 


h 


525.4 


523.98 


524.47 


524.60 


2 


21 


fi 


539.8 


540.92 


541.01 


541.27 


2 


22 



However, the apparent relative simplicity of a white dwarf 
means that a realistic model can be specified by only a few 
parameters and the potential richness of the g-mode pul- 
sation spectrum provides the possibility of many revealed 
modes, and hence many constraints to restrict the inversion 
process. 

Following a WET run on the object GD358 
l|Winget et al.lll994h th at yielded 154 h of near continuous 
time series photometry, iBradlev fc Wingetl (|l994l ) were able 
to use the detected pulsation modes to identify a preferred 
model for this star. They were aided in this task by the 
fact that the 11 detected normal modes clearly formed a 
sequence of increasing radial order (n) that were consistent 
with an i = 1 assignment, due to the evidence of a triplet 
structure for each mode. Their asteroseismic analysis en- 
abled a determination of such model parameters as: total 
mass, surface helium layer mass and effective temperature. 
Luminosity and a bolometric correction were also estimated, 
which then led to an asteroseismic distance determination. 

A further advance in these modelling procedures has 
been made more recently by the use of a genetic algorithm 
to efficiently explore the multidimensional parameter space 
in a global search for t he optimum model or models tha t 
fit the pulsatio n data ijMetcalfe. Nather fc Wingetl lioobi ). 
iMetcalfg (|2003l ) (and references therein) reports an anal- 
ysis of GD 358 data that is sensitive to the presumed white 
dwarf's 12 C/ 16 core composition, and thereby yields an 
indirect measurement of the 12 C(a,7) le O reaction cross sec- 
tion at astrophysically relevant energies (as the C/O core 
composition is determined by the competition between this 
reaction and the C-fo rming triple a reaction). Furthermore, 
iMetcalfe et ail (|2005l ) have successfully exploited the global 
optimization method to analyse pulsation data obtained 
from dual-site observations (and previous work) of another 
DBV star CBS 114, and clarify a number of features of its 
structure. 

Note that there has been some controversy con- 
cerning these asteroseismic successes, which is related 
to the fact that both composition variation in the 
core and in the envelope below the helium atmo- 
sphere pr oduce similar non uniformities in the period 
spectrum dFontaine fc Br assard 2002; Bras sard fc Fontaine! 
120031 : iMontgomerv. Metcalfe fc WingetJl2003l l~ 

Given that the analysis in the previous subsection did 
not uncover any clear trends in the pulsation mode data 
(except for perhaps the rotational and magnetic splitting 
we have proposed), one would hope that use of a global 
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fitting procedure could yield some definitive results. Adopt- 
ing this mode splitting model reduces the 11 identified nor- 
mal modes to only 8 independent modes with different n 
and/or I values. This is comparable to the number of modes 
available for the analyses of both GD358 and CBS 114 dis- 
cussed above. But, in contrast to these other two pulsators 
where additional constraints are evident (I = 1 assignments 
with consecutive radial index values), the EC 20058 pulsa- 
tion spectrum offers none of these clues. The revealed modes 
are most likely a mixtures of i — 1 and 1 = 2. 

Using the 8- mode data set, we applied a modified ver- 
sion of t he global model-fittin g procedure originally de- 
scribed bv lMetcalfe et al.1 (|2003h . This version of the code in- 
corporates the OPAL radiative opacities Jlglesias fc Rogers 
1996) rather than the older LAO data (jHuebner et al. 
19771 ). which are known to produce systematic error s 
in the derived temperatures l|Fontaine fc Brassard! 1 1994 ). 
The fitting procedure us e s a p arallel genetic algorithm 
l|Metcalfe fc Charbonneaul 120031 ') to minimize the root- 
mean-square residuals between the observed and calculated 
periods (op) for models with effective temperatures (T e ff) 
between 20 000 and 30 000 K, and stellar masses (K'h) be- 
tween 0.45 and 0.70 Mq. We restricted the mass range more 
than in earlier applications to avoid a family of models with 
high masses, which contain such a high density of £ — 2 
modes that they can match essentially any set of observed 
periods. We allowed the base of the uniform He/C enve- 
lope to be located at an outer mass fraction log(M env /M«) 
between —2.0 and —4.0. The base of the pure He surface 
layer could assume values of \og(Mn a /M±) between —5.0 
and -7.0. 

Since we have almost no information about the spheri- 
cal degree of the modes, we assumed only £ = 1 and I = 2 
modes were observable and calculated all of the periods in 
the range 150—600 s for each model; we then selected the 
closest model period for each observed mode. This proce- 
dure tends to bias the identification in favor of £ = 2 modes 
since t here are always more of them for a given model. Fol- 
lowing [MetcaifeetaD (|2004h . we required an £ = 2 mode to 
be closer to the observed period by a factor (A^ = 2/A^ = i) for 
selection as a better match. In effect, we optimized the mode 
identification internally for each model evaluation, while the 
genetic algorithm optimized the values of the other four pa- 
rameters. 

To quantify the effect of our ignorance of the core com- 
position, we repeated this fitting procedure using three dif- 
ferent types of cores: pure C, a uniform 50:50 mixture of 
C/O, and pure O. The results of these three fits are shown 
in Table [4] and the corresponding model periods and mode 
identifications are shown in Table [5] With the exception of 
the period at 350.6 s, the mode identifications are the same 
for all three fits. The fit using a pure C core identified the 
350.6 s period as (£ = 1, n — 7) while the other two fits both 
preferred an identification of (£ = 2,n = 13). The pure C 
model also includes this £ = 2 mode, and although it was 
closer to the observed period than the £ = 1 mode, it was not 
close enough to overcome the requirement that it be closer 
by the factor (Nt=2/Nt=i). Thus, the mode identification 
appears to be reasonably robust. 

All three core compositions yield very good fits to the 
observed periods, with op ~ 2 seconds in all cases. This 
is slightly better than the quality of the 4-parameter fits 



achieved for th e DBV white dwarfs CD 358 and CBS 114 by 
iMetcalfe et ail (|2005l) using the same code. The main effect 
of the different core compositions is to modify the optimal 
value of T e g, and to make slight adjustments to the loca- 
tions of the two near-surface composition gradients. The 
average values of the mass and effective temperature are 
in good agreement with the sp ectroscopically determined 
values (jBeauchamp et all Il999l ). the total envelope mass 
is within the range expected from stellar evolution theory 
(ID'Antona fc Mazzitehll 19791 ), and the thickness of the pure 
He surface layer, when compared to similar measures for 
CD 358 and CBS 114, is consist ent with the expectations of 
diffusion theory (|Metcalfell20"07l ). 



5.4 Review of the global fitting 

At this point it is prudent to pause and critically examine 
where our analysis has led us. From a total of 11 detected 
pulsation modes, we assumed that 8 of them had indepen- 
dent n and £ values, and then our modelling has identified 
unique index values for all but one of these modes (Table [5]). 
An obvious question is how plausible and unique is our 'best' 
model? 

One interesting result, which was foreshadowed in our 
previous general discussion, is the assignment of different i 
values (1 and 2) to the dominant 281s (fe) and 257s (fg) 
equal (observed) amplitude modes. Taking into account the 
expected larger geometrical flux cancellation from an £ = 2 
mode, this requires a higher physical amplitude for the 257 s 
mode to compensate for this geometrical effect. In isolation, 
a more plausible assumption might be that both modes have 
the same I — 1 value, and the small period spacing between 
them is explained by an unusually high mass (~ 0.8 Mq) for 
the star. But, we do have external constraints: model fits to 
EC 20058's optical spectrum are not consistent with suc h a 
high mass (|Beauchamp et al.l 1 19991 : ISullivan et alj 120071) . It 
is also interesting that a high mass can be ruled out largely 
through asteroseismic constraints, as discussed in the next 
paragraph. 

As a further check on our conclusions, we undertook a 
relatively exhaustive exploration of model parameter space 
using a grid of pre-computed models. Unconstrained fitting 
to the pulsation spectrum using criteria similar to that dis- 
cussed previously yielded two competing best fit models: 
the 'low-mass' (~ 0.55 Mq) ones discussed previously, and a 
'high-mass' family with M ~ 0.8 - 0.9 Mq. However, when 
the 281s (fe) mode was constrained to be £ — 1, then the 
high-mass family did not remain competitive, and hence the 
best model assigned £ — 2 to the 257 s mode and predicted 
a mass of 0.55 Mq, confirming the above. Assigning £ = 1 
to the 281 s mode is not totally arbitrary and is consistent 
with the rotational and magnetic mode splitting arguments 
developed in the previous section. 

Thus, two independent constraints added to the astero- 
seismic fitting resulted in the same conclusion. We conclude 
our modelling discussion at this point, and add the com- 
ment that a more exhaustive modelling effort is currently 
underway. 
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6 SUMMARY AND FUTURE WORK 

The work presented here has significantly extended our 
knowledge of the helium atmosphere pulsating white dwarf, 
EC 20058—5 234 beyond t hat re vealed by the discovery ob- 
servations of iKoen et al.l (1 19951 ). Analysis of nearly 135 h of 
photometric time series data obtained by four southern tele- 
scopes during the 8-day 1997 WET run XCOV15, along with 
42 h of single site data obtained during a week of observing 
at Mt John Observatory in June 2004 has more than dou- 
bled the number of detected frequencies from the 8 to 18 
(Table 

Following a detailed simulation procedure that we un- 
dertook, we are confident at the 1 in 1000 false alarm proba- 
bility level that the power in the detected modes (the lowest 
0.6 mma frequencies, in particular) is not simply the result 
of random noise fluctuations. 

Although we were able to identify signal power in our 
WET data at 18 different frequencies that are above our 
adopted 0.54 mma significance level (and one possibility just 
below: 'a'), only 11 of these could reasonably be viewed as 
pulsation modes. The 7 other periods are readily charac- 
terised as combination frequencies resulting from the non- 
linear processes that translate mechanical movement of the 
stellar material into observable flux changes. 

Our initial investigation of the detected pulsation mode 
structure revealed no clear systematic trends, except for per- 
haps a rotational splitting value of 70 (or 75) /iHz for two 
of the modes. Interestingly, the modes suspected of show- 
ing rotational splitting were independently identified as the 
only 1=1 modes from our objective fitting method, which 
could explain why they are the only modes that show this 
behavior. We (tentatively) interpret this as evidence that 
the white dwarf is rotating with a period of about 2 hours. 

A global fitting analysis of the period st ructure, follow- 
ing pr ocedures similar to those described in iMetcalfe et all 
(2000), resulted in 3 models with different core compositions 
that all yielded very good fits to the 8 assumed indepen- 
dent observed periods (see Tables [4] and [5} . All models had 
similar effective temperatures (~ 28 200 K) and masses (0.55 
Mq) that were in go od agreement with the spectroscopically 
deter mined values (IBeauchanrp et all Il999l ; ISullivan et al.l 
120071 ). This analysis provided no insights on the core com- 
position, as the fits were largely insensitive to the 3 different 
chosen C/O ratios. As expected from the earlier discussions, 
the global fitting assigned different mode I values to the two 
dominant pulsation frequencies: fe, I = 1 and fg, I = 2. 

Although the number of known DBV objects is rel- 
atively small (17), and is totally dwarfed by the current 
known number of cooler hydrogen atmosphere (DAV) pul- 
sators (now exceeding 140), it is still tempting to draw 
some conclusions about the nature and ext ent of the DB in- 
stabili ty strip a long the lines discusse d by IMukadam et alj 
l|2004l ) (see also IMukadam et al.1120051 ). Even using just the 
two stars GD 358 and EC 20058, we can infer that the 
DBVs mimick the characteristics exhibited by the DAVs. 
GD 358 has a rela tively low effective temperature (24 900 K, 
iBeauchamp et ahl ) and exhibits both large amplitude pulsa- 
tions and unstable behaviour (modes come and go), while 
EC 20058 has a significantly higher effective temperature 
(28 400 K) and displays low amplitude very stable behaviour. 
Similar behaviour is found in the ensemble of DAVs, and 



some authors have divi ded them into hot (hD AV) and cold 
(cDAV) pulsators (e.g. IMukadam et alJliooi) . The hDAVs 
feature only a few detectable pulsation modes that have low 
amplitudes and are relatively stable, while the cDAVs ex- 
hibit many variable modes. This is surely telling us some- 
thing about the onset and cessation of pulsation as the white 
dwarf cools through the temperature dependent instability 
strip. 

Certainly compared with GD 358, EC 20058 is a stable 
pulsator in terms of both mode amplitude and period. The 
two lower panels in Fig. [5] in particular, demonstrate that 
the mode amplitude stability might be subject to some qual- 
ifaction, as the mode f? has an amplitude of 1.2 mma in the 
WET data, but there is no evidence of it in the MJ04 data. 
An alternative possibility is that beating between unresolved 
modes explains the low-level amplitude variability. However, 
the dominant modes are always present in the photometry 
with about the same amplitudes. 

Both data sets presented here indicate that EC 20058 
has very stable pulsation periods. Thi s conclusion is also 
supp o rted by other observa tional work l|Sullivan fe Sullivan! 
l2000l : ISullivanl 120031 , 120051 ) . This means that we can at- 
tempt to measure evolutionary cooling effects over human 
timescales. The key point is that the DBV white dwarf 
model periods incr ease with decreasing temperature (e.g. 
iBradlev et al.| [l994). so careful monitoring of a stable mode 
period or periods over a multi-year interval should in princi- 
ple lead to a detectable period change and therefore a direct 
measure of the stellar evolutionary timescale. 

As one might expect, the models for both the cooler 
DAV white dwarfs and the much hotter DOV pulsating 
pre-white dwarfs also predict period increases with decreas- 
ing temperature. More than 10 years of time-series pho- 
tometry (including several WET runs) on the DOV ob- 
j ect PG 1159—035 produced a measured period increase 
(|Costa et al.lll999l ) , largely because the predicted rate is sev- 
eral orders of magnitude higher than for the DBVs. Note, 
that as a direct testimony on the difficulties of these mea- 
surement s, a 'false alarm' re s ult w as published at an ear- 



lier date llWinget et al.l|l985l . Il99 ll) . Kepler and collabora- 
tors (|Kepler et al.l 1005b ) have been monitoring the highly 



stable DAV G117— B15A over three decades and have fi- 
nally detected a definite pe riod change. As recently pointed 
out bv lWinget et al.1 (|2004l ). with sufficient quality data one 
might expect to measure a period change for a DBV over a 
timescale as short as ~ 5 years. This is the ai m of an ongoin g 
observational campaign featuring EC 20058 <|Sullivarl 2005). 

The cooling mechanism for both the very hot pre-white 
dwarfs and the hot white dwarfs is dominated by neutrino 
emission from the high-temperature degenerate cores (e.g 



emissi on trom the mgh-tcmpcraturc degenerate core s (e.g. 
Winget et aI1ll983l: lO'Brien fe Kawalerll200bl IWinget et at] 



200 



Models show that this core neutrino flux exceeds 
the surface photon flux down to an effective temperature 
of about 25 00 K, dependi n g on the mass of the star. 

In fact, IWinget et al. (2004) demonstrate that a DB 
model with a typical mass of 0.6 Mq with an effective tem- 
perature of 28 400 K (the EC 20058 estimated value), has a 
neutrino cooling flux approximately four times that of the 
photon flux. These facts suggest a very interesting experi- 
ment: measure a period change and compare this rate with 
two values obtained from models that include and exclude 
the neutrino mechanism. Then, provided we are confident 
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that the star's mass and effective temperature are in the 
right region, we have the exciting possibility that we are 
using a pulsating star to directly test a key low-energy pre- 
diction of electroweak theory. 

This is the current promise for EC 20058 as at these ef- 
fective temperatures the evolutionary cooling occurs at ap- 
proximately constant radius, in contrast to the much hot - 
ter pre-white dwarfs l|Winget. Hansen fe Van Horn! 1 19831 ). 
In the latter case, both continuing gravitational contraction 
and evolutionary cooling have an impact on the pulsation 
periods so any comparisons with theoretical predictions are 
less certain. 

Given the potential importance of this star, we have 
recently undertaken both time-series photometry and spec- 
troscopy using one of the Magellan 6.5 m telescopes 
l)Sullivan et al.l l2007). in a ddition to seaso nal monitoring us- 
ing Mt John Observatory (|Sullivanll2003l ). 
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